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We investigate the pumped charge and spin at zero bias by a modulation of two control parameters 
using the full counting statistics with quantum master equation approach. First we study higher 
order effects of the pumping frequency in general Markov systems and show in this limit the equiv¬ 
alence between our approach and the real-time diagrammatic approach. An adiabatic modulation 
of the control parameters induces the Berry-Sinitsyn-Nemenman (BSN) phase. We show that the 
origin of the BSN phase is a nonadiabatic effect. The pumped charge (spin) is given by a summation 
of (i) a time integral of the instantaneous steady charge (spin) current and (ii) a geometric surface 
integral of the BSN curvature, which results from the BSN phase. In quantum dots (QDs) weakly 
coupled to two leads, we show that (i) is usually dominant if the thermodynamic parameters are 
modulated although it is zero if the thermodynamic parameters are fixed to zero bias. To observe 
the spin effects, we consider collinear magnetic fields, which relate to spins through the Zeeman 
effect, with different amplitudes applying to the QDs and the leads. For interacting one level QD, 
we calculate analytically the pumped charge and spin by modulating the magnetic fields and the 
coupling strengths to the leads in the noninteracting and strong interacting limits. We show that 
the difference between these two limits appears through the instantaneous averages of the numbers 
of the electron with up and down spin in the QD. For the quantum pump by the modulation of the 
magnetic fields of the QD and one lead, the energy dependences of linewidth functions, which are 
usually neglected, are essential. 

PACS numbers: 05.60.Gg, 72.25.-b, 73.63.Kv, 05.70.Ln 


I. INTRODUCTION 

In a mesoscopic system, even at zero bias, a charge or 
spin current is induced by a slow modulation of two or 
more control parametersip— . This phenomenon, called 
quantum adiabatic pump, is theoretically interesting be¬ 
cause its origins are quantum effects and non-equilibrium 
effects. The quantum adiabatic pump is also expected to 
be applied to the single electron transfer devices and the 
current standar d^^d^ . 

The adiabatically pumped quantity is described by 
a geometric expression in the control parameter space, 
although the pumped quantity coming from second or 
more higher order of the pumping frequency is not ge¬ 
ometric. In noninteracting systems, the quantum adia¬ 
batic pump had extensively been studied by the Brouwer 
formulaM^— , which describes the pumped charge by the 
scattering matrix. On the other hand, it is difficult to 
calculate the scattering matrix in the interacting sys¬ 
tems. In the interacting system, the Brouwer formula 
had only been applied in mean field treatment a^^d^ or in 
the Toulouse limit^^. 

Recently, the quantum pump in interacting systems 
have been actively researched. There are three theo¬ 
retical approaches. The first is the Green’s function 
approach to pumping2^“— . The second is the real¬ 
time diagrammatic approach— (RT approach) which 
uses the generalized master equation (GME) that is 
equivalen t^^d^ to the quantum master equation (QME) 
derived using the Nakajima-Zwanzig projection operator 
technique^. Particularly, Ref. [s^ derived a geometric 
expression similar to the Brouwer formula and the Berry- 
Sinitsyn-Nemenman (BSN) vector explained later. The 


third is the full counting statistical^— (ECS) with quan¬ 
tum master equation (FCS-QME, which is also called 
generalized quantum master equation^) approach pro¬ 
posed in Ref. 

The adiabatic modulation of the control parameters 
induces a Berry-phase-like^ quantity called BSN phase 
in the FCS-QME with the Markov approximation. Sinit¬ 
syn and Nemenman^ studied the adiabatically pumped 
charge using the ECS and had shown that it is charac¬ 
terized by the BSN vector, which results from the BSN 
phase. The BSN vector was applied to the spin bo¬ 
son system^ and a connection was made to the excess 
entropy productioni2ii^. The FCS-QME approach can 
treat the Coulomb interaction, which can not be treated 
in the Brouwer formula. The derived formula of the BSN 
vector depends on the approximations used for the QME. 
The Born-Markov approximation with or without the ro¬ 
tating wave approximation^(RWA) is frequently used. 
The QME in the Born-Markov approximation without 
RWA sometimes violates the non-negativity of the sys¬ 
tem reduced density operator—. The QME of the RWA 
or the coarse-graining approximatio n^°i^^ (CCA) is the 
Lindblad type which guarantees the non-negativity^. 

Some recent papers^^iMii^ showed that the Coulomb 
interaction induces the quantum pump. In Refs. 
it was shown that in a one level interacting quantum dot 
(QD) weakly coupled to two leads, the pumped charge 
(also spin in Ref. [s^) induced by an adiabatic modula¬ 
tion of the energy level of the QD and the bias between 
the two leads vanishes in the noninteracting limit. In par¬ 
ticular, Yuge et ai— studied the pumped charge coming 
from the BSN curvatures by adiabatic modulation of the 
thermodynamic parameters (the chemical potentials and 
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the temperatures) in spinless QDs weakly coupled to two 
spinless leads and showed that the BSN curvatures are 
zero in noninteracting QDs although they are nonzero for 
finite interaction. 

In this paper, we first generalize the FCS-QME ap¬ 
proach to multicounting field to calculate spin current 
f ljll All . We then study the nonadiabatic effects in gen¬ 
eral Markov systems and clarified the relations between 
the FCS-QME approach and the RT approach^ in IIIIBI 
Additionally we show that the origin of the BSN phase 
is a nonadiabatic effect. Next we explain the model 
to be considered (flni). We consider QDs weakly cou¬ 
pled to two leads (L and R). To observe the spin ef¬ 
fects, we consider collinear magnetic fields, which relate 
to spins through the Zeeman effect, with different am¬ 
plitudes applying to the QDs (i?s) and the leads {Bl 
and Bji). The dynamic parameters {Bs, B^f^, and the 
coupling strengths to the leads) are control parameters. 
We use the RWA defined as a long coarse-graining time 
limit of the CGA to the FCS-QME. In llIVI and we 
consider noninteracting and interacting QDs respectively. 
First, we show (in l lIVBl and WBI) the time integral of 
an instantaneous steady current is usually dominant if 
the thermodynamic parameters (the chemical potentials 
and the temperatures of leads) are modulated (as consid¬ 
ered in Refs. dHlliisisi). Next in a one level QD 
with the Coulomb interaction U, we analytically calcu¬ 
late the BSN curvatures of spin and charge induced by 
the dynamic parameters in the noninteracting fl llV CD 
and strong f W Cl) interacting limit {U —>■ 0,oo). The 
difference between the results for 17 = 0 and U = oo 
appears through the instantaneous averages of the num¬ 
bers of the electron with up and down spin in the QD. 
For the quantum pump by the adiabatic modulation of 
{Bl,Bs), the energy dependences of linewidth functions, 
which are usually neglected, are essential. In HV Cl we 
show and discuss the contour plots of BSN curvatures 
evaluated numerically. Finally, we summarize this pa¬ 
per with discussions l UVIj) . In Appendix [Aj the Liouville 
space^2i^ and the matrix representation of the Liouvil- 
lian are explained. In Appendix [B] we derive the FCS- 
QME of the CGA and discuss the difference between the 
CGA and the RWA. In Appendix [Cl we derive Eq. (la¬ 
in Appendix |D1 we discuss the validity of the adiabatic 
expansion in IIIIBI In Appendix [El we discuss the deriva¬ 
tion of Eq. (IHTl) . 


II. FCS-QME 

In this section, we consider general Markov systems 
weakly coupled to noninteracting (fermionic or bosonic) 
baths. The model we use to do a concrete calculation is 
explained at HIIII 

In llll A1 we explain the FGS-QME method using the 
Liouville space^2i^(Appendix|A|. This method is a gen¬ 
eralization to the multicounting field of Ref. [i^. In 
1 |II Bl we study nonadiabatic effect, and show the equiv¬ 


alence to the method of Ref. . 

A. Derivation of FCS-QME 

Consider a cyclic modulation of the control parame¬ 
ters with a period r. At t = 0 and < = r, we per¬ 
form projection measurements of /ith time-independent 
observables {O^} indexed by p, of baths which com¬ 
mute with each other. Ao^ = o'^'^ — ofP denotes 
the difference of the outcomes at t = r and 

the outcomes at t = 0. The Fourier transform 

of the joint probability density distribution Pr({Ao;_j}), 

generating function. Here, Xi^ counting fields for O^. 
ZriiXt^}) is given by Zr{{xp}) = Tvtot[Ptot{t = ^)] using 
an operator of the total system Here, y de¬ 

notes the set of the counting fields {y^}. We defined 

p^{t) TrB[p((jj(t)] where Tr^ denotes a trace over 
baths’ degrees of freedom. p^{t) provides the generat¬ 
ing function Zr{{Xfi}) = Ti's[p^(^ = ''")]■ In Appendix iBl 
we derive the full counting statistics with quantum mas¬ 
ter equation (FCS-QME) [i.e., the equation of motion of 
p^(t)] from the equation of motion of p^^^{t). In this pa¬ 
per, we set h = 1. We suppose ptot(O) = p(0) C) Psioio) 
and ptot{t) ~ pit) ® psiat) (0 < t < r) where ptot(i) 
is the total system state, p{t) = TrB[ptot(Q] is the sys¬ 
tem reduced density operator, psioi-t) is a tensor product 
of the grand canonical or canonical distributions of the 
baths, and at is the value of the set of the control pa¬ 
rameters at time t. If the baths are electric leads, psiott) 
is given by Eq. (iBll . The FCS-QMEi2<43 jg 

^^=KX(^at)pHt), ( 1 ) 

and the initial condition is p^(0) = p(0). Here K^{at) 
is the Liouvillian modified by y. The Liouvillian de¬ 
pends on used approximations, for instance, the Born- 
Markov approximation without or within RWA^S and 
the CG A^°i^^ . After 91111 we choose the Born-Markov 
approximation within RWA; however, in this section we 
assume only Markov property (i.e., just depends on 
at). At y = 0, the FCS-QME becomes the quantum 
master equation (QME) 

^=kiat)pit). ( 2 ) 

K{at) equals K^iat) at y = 0. In the following, a symbol 
X without y denotes A^|^=o- 

In the Liouville space (Appendix |^, the left and right 
eigenvalue equations of the Liouvillian are 

k^ia)\pl{a))) = X^{a)\p^{a))), (3) 

ma)\k^ia) = XUa)ma)\. ( 4 ) 

The left eigenvectors l^ia) and the right eigenvectors 
Pmi^) (operators considered as elements of a vector 
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space) satisfy ((^^(Q;)|p))j(a))) = Snm- Here, a denotes 
arbitrary values of the set of the control parameters. The 
mode which has the eigenvalue with the maximum real 
part is assigned by the label n = 0 and is called the 
slowest mode. In the limit x —t 0, Aq (a) becomes 0 and 
((Zg(Q;)| becomes ((1|, i.e., lo{c() = 1. The conservation of 
the probability ^((l|p(t))) = {{l\K{at)\p{t))) = 0 leads 
((l|7ir(a) = 0. In addition, n = 0 mode right eigenvector, 
|po(ck))), determined by K{a)\pQ{a))) = 0 represents the 
instantaneous steady state; if the control parameters are 
fixed to a, the state p(t) converges to po{oi) at t —>■ oo. 
In general, the solution of the FCS-QME Eq. m is ex¬ 
panded as 


\pHt)))=Y,cme^-^*HpUat))), (5) 

n 


where A^{t) = fg ds A^(as). The coefficients c^(t) obey 
Eq. (IC8l) . 

The coefficients c^(t) are given by solving Eq. (IC8I) . 
The condition which makes {c))(t)}„/o negligible is dis¬ 
cussed in Appendix [C] and was studied carefully in Ref. 
[S^ for X = 0- In this section, we consider sufficiently 
slow modulation of the control parameters. The effects 
of fast modulation are considered in the next section. For 
the slowest mode, the second term of the right side of Eq. 
(ICSp exponentially damps as a function of time. The re¬ 
laxation time of the system (rs) is the order of F“^ where 
F is the typical value of the linewidth functions [defined 
as Fi in Eq. mi Assuming the cycle time r is much 
longer than rg, we obtain 


c^(^) = c^(0)exp 


dt miat)\-\p^{at))) 


,( 6 ) 


and 


~ 0 (n ^ 0). (7) 

In fact, = O(^) with w = 27r/T as explained 

in Appendix [C] We denote the condition which makes 
the above approximation appropriate as x-O'diabatic con¬ 
dition. Using the initial condition p^(0) = p(0), we 
obtain Cg(0) = ((Zg (ao)|p(0))). Substituting these equa¬ 
tions into Eq. we obtain^ 

~ {{lo{ao)\p{0)))e- dS 

( 8 ) 

and the cumulant generating function >S't(x) = 

In^r(x) = ln((l|p^(T))) ; 

SAx) = £dtXUat)-£da- 

+ ln((Zg^(ag)|p(0)))+ln((l|pJK))). (9) 

Here, we used Jg dt ((^g (aOIJ^I/Og (a*))) = 

ff^da^ ((^0 where C is the trajectory 


from ao to ar, a" are the nth component of the control 
parameters and the summation symbol ^ is omitted. 
Equation (0 is the same with Yuge et alM except for 


that X denotes a multicounting field. 
(Ao^)t- = are 


The averages 


(Ao^),- = J dt X^{at) - <j> da^ ((^g (a)| 

+((/g^K)|p(0))) + ((l|p^(ag))), 


5|Po(a))) 


da" 


( 10 ) 


where X^(a) Ix-o' Here, we used 

-§^da" {{lg{a)\^^^ = -((1|p(;K))) + ((l|p(;(«o))) 

because ((Zg(a)I = g| 7 r((l|Pg («)))• The integrand 

of the first time integral, Ag(at), are the instantaneous 
steady currents of at time t\ if the control parameters 
are fixed to a and the state is po{oi), the current of 
is Xg(a). The integrand of the second term of the right 
side of Eq. m is the Berry-Sinitsyn-Nemeninan (BSN) 
vector— 


MJ./' ^ // 7 P/ M^I/^o(c^))) /"I 1 ^ 

^n('a) = ((Zo(a)l ■ (11) 

The third and fourth terms of the right side of Eq. CHD 
cancel if the initial condition is the instantaneous steady 
state pg(ag). Because of a,- = ag, the second term of 
the right side of Eq. (nni can be described as a surface 
integral over the surface S enclosed by C using the Stokes 
theorem : 

(Ao^). = (Ao^)f-'i^ + (Ao^)f-^ (12) 

(Ao^)f-'iy= Tdt A^(at), (13) 

(Ao^)f = - f da^ A da- (14) 

^ s 

Here, A is the wedge product and the summation symbol 
Sn m omitted. BSN curvature F)((jj(a) is given by 


Knic^) = 


dAl^ia) dAi^{a) 


da" 


da" 


(15) 


Yuge et alM focus on only the second term of Eq. m 
subtracting the first term, and they did not evaluate 
In 9IVB1 we show that this contribution is 
usually dominant if the thermodynamic parameters are 
modulated although the steady currents Xg(at) are zero 
if the thermodynamic parameters are fixed to zero bias. 

From the cumulant generating function Eq. (j9|), we 
can calculate the second order cumulants ®(o^Oi,)r = 

d(ix%[^x ) lx~o higher cumulants. However, we fo¬ 
cus only on the first order cumulants (averages) in this 
paper. Up to the first order cumulants, we do not need 
the multicounting field. In fact, the multicounting field 
is helpful to understand the physical origin of each term 
of the (FCS-)QME. The counting fields appear only in 
A» B (A, B 1) type terms in the dissipator term H^* 
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[see Eq. (IA7I) : • is an arbitrary operator] as Eq. (jB7l) . 
For instance, the factor with the counting field for 
total number of electron with spin a of bath &, means 

tunneling process from the system to bath b with spin a. 
The factor means tunneling process from bath b 

with spin a to the system. 


B. Non-adiabatic effect and BSN vector 

In this section, we consider nonadiabatic effect [that 
comes from (n = 1, 2, • • •) with w = 27r/T], which 

had been researched recently^iiS^i^i^. If the modulation 
of the control parameters is not adiabatic, the difference 
between the state and the instantaneous steady state, 
p(t) - po{at) = is im¬ 

portant. In contrast to Appendix IQ and Ref. [l^, we do 
not treat {c„(t)}„^o explicitly; instead, we use the pseu¬ 
doinverse of the Liouvillian. The formal solution of the 
FCS-QME Eq. (P) is 


((l|pg (ao))); result from this relaxation. The contribu¬ 
tion of (Ao^)t- from (5p(0) = p“(0) = p(0) — Po(q^o) is 

(Ao^r = 

= ((/';(ao)|dp(0))) - {r,[a.)\5p{T))) 

+ ( 20 ) 


with |(5/9(t))) Texp Jq ds K{as) |i5/o(0))). The first 


term of the right side of Eq. (l20ll is (Ao^)^+'^. Because 
we can obtain ((Zq (a)|/9o(a))) + ((l|Po('^))) = 0 from the 
normalization ((Zg (a)|/9g (a))) = 1, (Ao^)^+^ is given by 
((^o(«o)[|p(0)))-|p(ao)))] = ((^^(ao)|(5p(0))). The second 
term of the right side of Eq. (EQl) is exponentially small 
since 5p{t) ^ The order of the third term is C>(^) 

with bj = 21 :It because = 0{lu) and the integral 

range is restricted up to l/T since Sp{t) ^ e“^‘. Hence 

{Ao,yf = {Ao,)l+'^ + Oif). 

The currents can also be written as 


|pX(t))) = Texp / ds (as) \p^{0))), (16) 


'-Jo 


where T denotes the time-ordering operation. Using this, 
we obtain the averages^ 


4(<) = ((l|IU^(aO|p(t))), (21) 

where W^(a) are the current operators defined by 

((l|IU^(a) = {{l\k>^{a), (22) 


(Ao^)t = 


diixkl 


{wm 


x=o 


= du {{l\kf^iau)\p{u))) + ((l|p'"(0))) 

Jo 

= f du {{l\kf^{au)\p{u))) = f du Ifj,{u){17) 
Jo Jo 


Here, we use ((l|A'(a) = 0 in the second line and 
|P^(0))) = 0 (derived from |p^(0))) = |p(0)))) in the third 
line. Moreover, using ((Zo(Q')I = ((1|, •^o(a) = 0 and Eq. 
O, we obtain^ 


{{l\k^{a) = A[;(a)((l| - {{lika)\k{a). (18) 


i.e., Trs[IT^(a)*] = Trs[Ar^(a)*] for any operator • [see, 
Eq. (IA14I) ]. Using Eq. (fTSl) . the instantaneous steady 
currents are given by 

((l|lU^(a)|po(a))) = K(a) = /^-'^^(a). (23) 

In the QDs weakly coupled to two leads, the electric cur¬ 
rent operator [i.e., lU^(a) corresponding to the electric 
current] coincides with Ref. in the Born-Markov ap¬ 
proximation without or within RWA. 

Applying the pseudoinverse 7Z(a) defined by 

7^(a)A'(a) = l-|po(a)))((l|, (24) 

to the QME Eq. m, we obtain 


Hence the currents are given by 

4(t) = ((llk'^(at)lpm 

= K(a,)-((l^(a,)lk(a,)lp(t)}} 


(19) 


|p“(t))) = 7Z(a,)^lpo(a,))) + 7^(a*)||p“(^))) 

OO J CO 

= E [nat)^Y\po{at))) ^ IP“^"Hi)))(25) 

n—1 n—1 


Substituting Eq. (1251) to Eq. eo, we finally reach 


The used approximation is only the Markov property of 
the Liouvillian, Eq. O- In Appendix IdJ we discuss the 
reasonable range of the power of ^ (nonadiabaticity). 
Substituting Eq. (fTOl) with p{t) Ki po{cxt) into Eq. (IT71) . 
we obtain Eq. cni) without the third and fourth terms. 
If p(0) ^ Poicto)-: the state relaxes to the instantaneous 
steady state in the relaxation time ts- The third and 
fourth terms of Eq. (HU]), (Ao^)^+^ = ((?(( (ao)|p(0))) -I- 


n—1 

with ((l|IU^(at)|p“("')(t))). Since = 

C>(a;) and TZ{at)=0{^), 

P“(”)(t) = 0(p”- 


(27) 
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In Appendix |D1 we discuss the reasonable range of n of 
^a(n)(^) show that with the larger nonadiabaticity 
(p), the reasonable range becomes wider. 

Let’s consider the relation between Eq. dT51) and Eq. 
(l26l) . In gil A[ we used y-adiabatic approximation Eq. 
JSl), which becomes |p(t))) ~ \po{at))) at y = 0. Substi¬ 
tuting it to Eq. (1^ . we obtain I^{t) « So, 

we cannot obtain nonadiabatic currents ■ 

However, from the y^ derivative of Eq. ([5]), we obtain 

(28) 

This is equivalent to Eq. for p(0) = po(<ao)- Equa¬ 
tion (l28l) suggests 

i;^^Ht) = -{{li;{a,)\-^^\po{a,))). (29) 

In fact, this is equivalent to = 

{{l\Wfj,{at)\p°-^^\t))), namely 

/“(i)(t) = {{l\WMna,)f^\po{a,))), (30) 

because of 

((l|IE4a)7^(a) = -{{li;{a)\ + c^(a)((l|, (31) 

which was shown by Sagawa, et al^ for a single count¬ 
ing field. Here, Cfj^{a) are constants shown in Eq. (lE7l) . 
We prove Eq. pTI) at Appendix [E] Equation ([191) 

Eq. (l26ll are identical because of Eq. (l23l) and Eq. 
(I5T]) . In other words, in the expansion of I^{t) obtained 
from a substitution of p{t) = po{at) + Yl^=i P°‘^'^\t) = 
J2‘^=o info Ed- (HU)! tho nth (n = 0,1, • • •) or¬ 

der nonadiabatic solution, p“(")(t), gives {n + l)th or¬ 
der nonadiabatic currents (i) because of Eq. (ISTl) . 

Hence the ECS-QME approach picks out one higher or¬ 
der nonadiabatic piece of information from the solution 
of the QME. 

Moreover, although the BSN phase [i.e., the argument 
of the exponential function of Eq. (|6|)] is derived under 
the y-adiabatic condition which makes Eqs. m and (ED 
appropriate, its origin is probably a nonadiabatic effect 
that comes from because Eq. (1^ shows that the BSN 
phase has the information of the nonadiabatic part of the 
QME [p“(t) = p{t) - po{at)]. 

It is important to recognize the relations between the 
FCS-QME approach and the RT approach^—. In the 
RT approach, pK.{t) = (K|p(t)|«;) are governed by the gen¬ 
eralized master equation (GME) 


d 

dt 






dt’ 


(32) 


where |k) are the energy eigenstates of the sys¬ 
tem Hamiltonian. The kernel HQ,;(t,t') can include 
the higher order contribution of the tunneling in¬ 
teraction between baths and the system. In the 


GME, p,{t') is given by p,(t) + 

Moreover, Wnnitjt’) and Pn{t) are expanded as 

and p^(t) = 

Er=o ^ P^vU) 

are of the order of w"W. In particular, it]t — t') = 

(at; t — t') is the kernel where the control param¬ 
eters are fixed to at- Up to the second order of the tun¬ 
neling interaction (in the following we consider this level 
of approximation), we obtain^i*^ 

0 = ^iLW(at)p(°)(at), (33) 

r) 

= (34) 

V 

for n = 0,1, • • •, with 

iL(°)(at)=r dt’ (35) 

J —OO 

which is the instantaneous Liouvillian corresponding to 
our K{at). Equation (1551) is just the definition of the 
instantaneous steady state pl^\at) = p|^°p^(t), which sat¬ 
isfies = 1- Additionally, P^^(j){t) for n > 1 

satisfies E«pLo) {t) =0. The conservation of the prob¬ 
ability leads to (at) = 0, which corresponds to 

our ((I|A'(at) = 0. The charge or spin current /^(t) is 
given by2^i2^ 

= (36) 

K,r) 

corresponding to our Eq. EID- wL^(at) is the instan¬ 
taneous current matrix of in the present approxima¬ 
tion, which corresponds to our W^(at) and is linear in T. 

Substituting pr,{t) Ri Eq. dMD! we 

obtain 

OO 

ut) = ^^^(^) = E^!((](«*)4?h)(^)(37) 

n—O K,ri 

Equation (p4l) for n = 0 leads to^ 

= (33) 

K 

Here, i?^„(at) is the pseudoinverse of (at) corre¬ 
sponding to our TZ(at) and it is given by^ 

i?^„(at) = (K-\., Kr,. = (39) 

Substituting Eq. ([38)) into Eq. (|57|) . we obtaiit^ 

4;‘>«) = i:Aka.)!^|)^. (40) 

K 

7’L'''(at) = E^C^'(«i)-^')34(«‘)- (41) 
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A similar method has been used in Ref. [s^l. (at) and 
Eq. (HOI) respectively correspond to our ((l|lE^(a)7?.(a) 
and Eq. dsnu- Moreover, Eq. (1341) for arbitrary n leads 
to 

PTii-n-l)W ~ ^ , (42) 

K 

which corresponds to our Eq. ([25|l . Because of these re¬ 
lations, the RT approach is equivalent to the FCS-QME 
approach in the calculation up to the second order of the 
tunneling interaction. Additionally, we discuss correc¬ 
tions due to the nonadiabatic effect of the FCS-QME in 
Appendix [P] The first equation of Eq. (ID6I) is consistent 
with p|j^(o)(t) = 0{ujtb), which can be derived from Ref. 
M- Here, tb is the relaxation time of the baths. 

In this section, we proved the equivalence between Eq. 
(HH) and Eq. (E51) using a key relation Eq. (1511) and 
showed the origin of the BSN phase is a nonadiabatic 
effect, and connected the FCS-QME approach and the 
RT approacbSi. These are among the most important 
results of this paper. 

III. MODEL 

We consider quantum dots (QDs) (denoted by a sym¬ 
bol S) weakly coupled to two leads. The total Hamilto¬ 
nian is Htot{t) = Hs{t) + Y.b=L,RiHb{t) + Hsb{t)]. Here, 
Hs{t) is the system (QDs) Hamiltonian, i?f,(t) is the 
Hamiltonian of lead b = L,R, and Hsb{t) is the tun¬ 
neling interaction Hamiltonian between S and lead b. To 
observe the spin effects, we snppose that the leads and 
the system are applied to collinear magnetic fields with 
different amplitudes, which relate to spins through the 
Zeeman effect. The leads are noninteracting: 

^b{t) = ^{sbk + crgbBb{t))cl^^Cbka- (43) 

k,a 

Here, a =t, 4.= ±1 is spin label, gb = where is 

the g-factor of lead b, fiB is the Bohr magneton and i?b(t) 
is the strength of the magnetic field of lead b. cl/.^(cbkcr) 
is the creation (annihilation) operator of an electron with 
spin (T and momentum k in lead b. The system Hamilto¬ 
nian is 

^ns,ms' (Bsit)) T .Hooulombi (44) 

n^m,s,s' 

where is the creation operator of an electron with or¬ 
bital n and spin s. £ns,ms'{Bs{i)) means the energy of the 
electron for n = m, s = s' and the tunneling amplitude 
between orbitals for (n, s) ^ {m, s') which depends on the 
magnetic field of the system. i?couiomb denotes Coulomb 
interaction. The tunneling interaction Hamiltonian is 

B^SbijP) — ^ ^ \/ Af, i^Vbkiy^nsbb\isb^bka -(-h.C., (45) 

k,(T,n,s 


where Ah(t) is a dimensionless parameter, and Vbka,ns is 
the tunneling amplitude. 

We assume Rg, Bb/r and A ^,/r are control parameters 
(denoted a' = {Rs, Rl//j, A^/zj} and are called the dy¬ 
namic parameters). The thermodynamic parameters (the 
chemical potentials and inverse temperatures of leads, 
{/ife} and {fib}) are also considered as control parameters 
in mVBI and HVBl We denote a" = {/3h,/ib}b=L,_R and 
a = a' + a". Yuge et alM chose the set of control pa¬ 
rameters as only a". However we are interested in a' for 
the reason explained in llIVBl 

We choose the measured observables {O^} = 

^ J2kbika<^bka- The pumped 
charge (spin) of lead b is given by (ANt^) ± (AiV^j,). 
{ANba) are calculated by Eq. (IT^ . In fact, what 
we call the pumped charge, (AWt) + (AA^j,), is the 
pumped electron number (actual pumped charge is given 
by —e[(AAf,.t') + (AAb 4 ,)], where e (> 0) is the elementary 
charge). 

In llIV Cl and (|V] we consider a one level system 

Hs{t)= ^ LOsiBs{t))alas+ Ua\.afalai, (46) 
s=t,4- 

as a special model of Eq. (|4^ . Here, s =1,4-= ±1, 
uJs{Bs) = Wo -I- sgsBs with wo the electron energy at 
Rs = 0, and gs = ^where g^ is the g factor of the 
QD. 

In the following, we apply the FCS-QME with rotating 
wave approximation (RWA) explained in Appendix |B1 


IV. NON-INTERACTING SYSTEM 

In this section, we consider a noninteracting system 
(Rcouiomb = 0). The system Hamiltonian Eq. (l44ll can 
be diagonalized 

2N 

Hs = J2^^blb,^ ( 47 ) 

i=l 

by a unitary transform a„s = Uns,ibi- The tunnel¬ 

ing interaction Hamiltonian Eq. (I45|) is 

Hsb = ^ Wbka,ib\cbko + h.C., (48) 

k.(T,i 

with. WfofcfTjZ — ^Zn,s "n/ ^b'f^bka-,nsU^g^j^. 

In llIV A[ the Liouvillian and its instantaneous steady 
state are explained. In llIVBl we consider the contribu¬ 
tion of Eq. (TT^ and show that this cannot be neglected 
in general if the chemical potentials and the temperatures 
are not fixed. In l lIV Cl we calculate the BSN curvatures 
for two combinations of modulated control parameters 
(Rl,Rs) and (Al,Rs). 
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A. Liouvillian 


The Liouvillian in the RWA is given by 

2N 

k^a) = Y^kf{a), (49) 

Kf{a)» = •] + nf (x, a) • +flf (a)*, (50) 

if {Qi} are not degenerated. Here, superoperators 
nf (x, Q:) and nf (a) operate to an arbitrary operator • 
as 


nf(x,a) 


E 


• b] - - 2‘^kblb^ • }> (51) 


nf («)• = i[ili{a)blbi,»], 

with 

k,(T 

xS{ebk + crgbBb - ul^), 


(52) 


1 


k,(T 


^bk H“ ^Qh^b 


(53) 

,(54) 


and n,(a) = ^ Eb • Here, /+(a;) = 

_|_ is i;]^e Fermi distribution function, 

fbi^) = 1 ~ /b'(^), Xb(7 is the counting field for Abo¬ 
und P denotes the Cauchy principal value. The ma¬ 
trix representation of kf{a) (see Appendix by the 
number states of blbi (|0)i and |l)i) is a 4 x 4 matrix 
which is block diagonalized to {|0)ii(0|, |l)ii(l|} space 
and {|0)ii(l|, |l)ii(0|} space. The {|0)ii(0|, |l)ii(l|} part 
is given by 




(55) 


with = Eb^M • {|0)«(1U1)**(0|} part does not 
relate to the instantaneous steady state of kf{a). The 
eigenvalue of the instantaneous steady state of kf{a) is 
given by 

^x ^*^(a)+ $,"(«) 

— r, 


+ jDf(a), (56) 


with D^(a) = [4>f -b j^/d - [4>f ’^4>f’^]. The 

corresponding left and right eigenvectors are |pfo('a))) = 
Cf(a)|00)). + E^(a)lllh and ((ll(a)l = ’.((00| + 


uf(a)j((ll| with Cf(a) = 
[Af +$+] 2 -i- 3 >"'>^<i.+-^ ’ and 




Af (a) = 


^ 4>+ - + 2/DfH 


2 $ 




(57) 


At XbCT = 0, Af (a) becomes Ei(a) = and Cf (a) 

becomes Ci (a) = 1 — Ei(a). 


B. Instantaneous steady currents 


The instantaneous steady current is given by 
^Steady (^) ^ noninteracting sys¬ 
tem, Ag(a) is Ei-^^o(“) and it leads to = 


Here, 

lated from Eq. (l56l) as 


Ix=0 


are calcu- 


^Steady 


{a) 


r£a,»rfly(/fl(c5t) - fL{oji)) 


(58) 


with Tbcr,* = 27rEfe \Wbka,i\'^b{ebk + crgbBb - cbi), Tb,i = 
Ecr^ba,i, and Ti = Ebl'bd- vanishes at zero 

bias {Pl = Pr, g-L = gn)- Let us consider the mod¬ 
ulation of only the thermodynamic parameters (a" = 
{gb, Pb}b=L,R) similar to Refs, . The fac¬ 
tor depending on a" of is - 

//3i,(t).A«i,(t)(w*)) with // 3 ,^(w) = + l]-i. Hence 


^ [// 3 K(t),AiK(t)(^») ~ // 3 i:,(t),Aii(i)(^ 0 ]( 59 ) 

is generally nonzero and is much lager than (AAbo-)^®'^'^^ 
because the period t is large for adiabatic pumps. Simi¬ 
larly, we can show that (AAf,^)®*®'^'^^ is generally nonzero 
for interacting system I WBI) . Reference considered 
special modulations of only thermodynamic parameters 
which satisfy = 0. In fact, the instanta¬ 

neous steady currents are always zero for arbitrary mod¬ 
ulations of only the dynamics parameters at zero bias. 

The pumped charge and spin due to the instantaneous 
steady currents (backgrounds) are generally nonzero even 
if the time averages of the bias are zero. References 
[mi chose V = gL — gn as one of the modulat¬ 
ing parameters and considered a pumping such that 
i fj dt V{t) = 0 and (AAh^)S‘«‘^dy ^ q. In such pump¬ 
ing, the (thermal or voltage) bias is effectively nonzero. 

Even if the backgrounds do not vanish, one can detect 
the BSN curvatures by subtracting the backgrounds by 
using zero-frequency measurements or by lock-in mea¬ 
surements. However, if one wants to apply the adia¬ 
batic pump to the current standar d^^d^ ^ the instanta¬ 
neous steady currents should be zero at all times because 
the backgrounds are sensitive to the velocity of the mod¬ 
ulation of the control parameters and its trajectory. In 
contrast, the pumped charge and spin due to the BSN 
curvatures are robust against the modulation of the ve¬ 
locity and the trajectory. Hence, if one wants to directly 
apply the BSN curvatures to, for instance, the current 
standard, one should fix the thermodynamic parameters 
at zero bias. 


^L<T,kR,i 


r* 


















C. BSN curvatures 


In the following, we consider one level system of which 
the Hamiltonian is Eq. dUl) at ?7 = 0. The instanta¬ 
neous steady state is given by I/Sq (o^))) = ®s=t.4-IPs o(“))) 
because the Liouvillian is described by a summation 
(Kx = Similarly, the corresponding left 

eigenvalue is given by ((/o(a)| = ®s=td((^?,o(“)l- The 
BSN vectors Eq. m are given by 


where 


S=t.t 


dEs{a) 

da" 


,.bcr 


{a') 


d{iXba) 


x=0 



(60) 


(61) 


with 


^ba,s{a') = 27rA6^ \Vbka,s\'^ 
k 

xS{£bk + crgbBb - LUo - sgsBs). (62) 


v^f{a') dose not depend on a". Equation (|60)l leads to 
an expression of the BSN curvatures 


t’Xh = E 


Ydv^f{a') dEs{a) 
I da"^ da" 


(m -n- n) . (63) 


We emphasize that Eq. (l63)) is consistent with the re¬ 
sults of Refs. which showed that the pumped 

charge (and also spin in Ref. [s^) vanishes at the non¬ 
interacting limit in these settings. The set of control 
parameters a was a" (for Ref. 0) and {wo, V = 
— Mfl} (for Refs. [^1^ 1. If a"* or a" is an ele¬ 
ment of a", F^^{a) is consistently zero. In Refs, 
the linewidth functions were energy-independent, namely 

rb(T,s(aO = i5(T,sr{,=constant. Hence ^ = 0 = 

° ^ ~ f* f'old consistently. 

To calculate F^^(a), we need to assume the energy 
dependences of r&o-^s- For the simplicity, we assume that 


r&cr.s = ^cr.slTb -I- r'f, • (sgsBs - crgbBb)] 

= Srj^sAbljb + 7b • {sgsBs - crgbBb)], (64) 


where TJ, are energy differential coefficients of linewidth 
functions at Bb = Bs = 0. Namely, we disregard spin 
flips induced by tunneling between the QD and the leads. 
Equation (IMl) is always appropriate when jT^gsHs — 
gbBb)\ <C Tb is satisfied. Additionally, we fix a" to zero 
bias {/3b = f3, gb = g), in which Eg/a) is given by Eg/a) = 
/(wq + sgsBs) with /(w) = -|- 1]“^. In this 

condition, (a™, a”) = {B^, Bs), (A^, Hg) components of 
the charge and spin BSN curvatures of lead L are 


^bI,Bs T ^Bt,Bs = -gsgL^Llfi^^o + 9sBs) ± f {ujq - gsBs)] 




^ tot 


+9sgLT's[f'{i^o + gsBs) T /'(wo - gsBs)] 

2Tr 


X (^'l(9sBs - gLBL)—^ + ^'r{9sBs - gRBR )-^— -Y 


Pkl,Bs ± ^At,Bs = 9s[f{uJo + gsBs) T /'(wo - gsBs)] T^Afl)2 


+5s[/'(wo + gsBs) ± /'(wo - 9sBs)]xl{9sBs - 9lBl) 


IrAr — "IlAr 
{IlAr + JrAr)^ 


(65) 


( 66 ) 


Here /'(w) = and Ttot = Tg -|- T^. The pumped 

charge (spin) induced by a slow cycle modulation of 
{a", Bs) {a" = Bl,Al) are given by 

(ANl^) ± (ANli) 

= - [ da-dBs ± fY\r,), (67) 

JS’' 

where S" are areas enclosed by the trajectories of 
{a",Bs). T ^a^,Bs = Bl,Al) are invariant 

under the transformation 7 b —>■ cy, —>■ cy( (for any c > 
0). Hence relevant quantities are y^/Ttot- The coupling 
strength Ttot itself is not important. ± Fg^ 


/ i" r \ 

are proportional to gsgL and E^)^ g^ ± PaI,Bs are pro¬ 
portional to gs- The first terms of the right side of Eq. 
(1651) and Eq. (l66)) are dominant terms. In the limit 
Yg —7> 0, EgYss T fYYbs and the second term of Eq. 
dSSl) vanish; however, the dominant term of Eq. ((Ml) 
remains. At ujo = g, f{uio + gsBs)-f'{uJo-gsBs) van¬ 
ish. Hence, at wq = the dominant terms of the spin 
BSN curvature of {Bg, Bs) pump and the charge BSN of 
{Ag, Bs) pump vanish. The contour plots of these BSN 
curvatures are shown in Figs. [Ija) and 1(b) and Figs. 
UKa) and 2(b). The details are explained in W Cl 

It is important to remark that {a"", a") = 
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(Si,(Ai, A/{) components of the charge and spin 
BSN curvatures are zero at zero bias because, in Eq. (1551) , 
Es{a) = /(wq + sgsBs) are independent of Bljji and 

^L/R- 

V. INTERACTING SYSTEM 

In this section, we study the interacting system Eq. 
dm. First, we explain the Liouvillian for 0 < C/ < oo 
f gyAl) . Next, the instantaneous steady charge and spin 
currents are calculated at U = oo f WBI) . In HVCl we 
confirm the consistency between our results and Ref. 
for 0 < 17 < oo. Finally, the BSN curvatures correspond¬ 
ing to Eq. (1551) and Eq. (1551) are calculated at C/ = oo 


and differences of the results between [7 = 0 and U = oo 
are discussed. 


A. Liouvillian 

We explain the Liouvillian for k^T > E, in which the 
Born-Markov approximation is appropriate. The matrix 
representation of the Liouvillian of RWA by the num¬ 
ber states = 0 i 1 a-re the numbers of an 

electron with spin s =t, 4 -) is a 16 x 16 matrix which 
is block diagonalized to the “diagonal” space (spanned 
by “off-diagonal” space 

(spanned by 

agonal” block is given by 


J 


K^{a) 


/-[$+ + $+] 

ci,+.x -[$-+</>+] 

0 

V 0 


with 

= 27rAt ^ \vbka,s\'^f^{uJo + sgsBs + U) 

k,a 

-h agbBb - wo - sgsBs - 17),(69) 

and = (j)'^’g^\u=o- The off-diagonal block is a 

(12 X 12)-diagonal matrix, which dose not relate to the 
instantaneous steady state. At 17 = 0, K^{a) becomes 
Ar^(a)® I 4 .-bit® at( a), where Ar 4 ^(a)(s =t, 4 ') are given 
by Eq. (1551) and 1^ are identity matrices. In the opposite 
limit f/ —>■ 00 , K^{a) reduces to 




/-[$+ + $+] 10000)) 

-<t>- 0 1010))(70) 

V ) 10101)) 


because the density of state of both leads vanishes at high 
energy {cj)f 0 ). 


B. Instantaneous steady currents 

In this section, we set U = 00 . The characteristic poly¬ 
nomial of A:x(oo) jg denoted as C' 3 (x, A) = det(A:x(°°) - 
^) = C'n(x)A" - A^. Because of (70(0) = 0, A = 0 
is one of the solutions at y = 0. Now we set Xba as 
infinitesimal and other counting fields are zero. Then, 
the eigenvalue corresponding to the instantaneous steady 
state is given by A = Ao(x) = iXba ■ -b 0{xl^). It 


0 

-[d>J+4] 

X,+ ’X 


0 \ 



-[•Pt +hV 


10000 )) 
11010 )) 
joioi)) ’ 
11111 )) 


( 68 ) 


leads to 0 = C'3(x, Ao(x)) = Ci{0)ixbal!^^^^ + iXb.C^'^ 
with Ck°' = I „, and we obtain 

^ ^(^Xbcr) IX=0’ 


rSteady 

^b<7 


Cl 


ba 


Clioy 


(71) 


with Cl (0) = — [$)!■ -b -b d>.j. ]. From Cq (x) = 

we have 


rSteady 

^Lg 


{a) 


Ss=t4 ^-B^L(y,s^R,s{fR,s - fL,s) 


,(72) 


where $ 1 ^ (s =t,-l') describes for s =t and for 
s = 4 ,. At zero bias, the instantaneous steady currents 
vanish. Similar to dlVBl are generally 

nonzero when a" is not fixed at zero bias. 


C. BSN curvatures 


The instantaneous steady state poi^a) and correspond¬ 
ing left eigenvector Ig (a) are written as po = po|00)(00| -b 
p^|10)(10| + p;|01)(01| + P 2 |ll)(ll| and ^ = |00)(00| + 
1^|10)(10| -b /^|01)(01| -b I 2 |11)(11|- The BSN vectors are 
given by 


A)('^(a)= ^ l^(a) 

c=td,2 


dpc{a) 

da^ 


( 73 ) 
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Charge( U=0) 


Spin([/=0) 


1.0 X 10- ' 
O.C X K)- ' 



-10 -5 0 5 

Charge( f/=Qo) 



-10 -5 0 5 

Spin(f/=Qo) 


-10 -5 0 5 

/'B^s/r 




-5 0 5 

//B^s/r 


?L1- 




LI 


]/( 


at C/ = 0, (b) the BSN curvature 

Z/'f '' ' ^ ^ ^ 


FIG. 1: (Color online)(a) BSN curvature of charge of (Bl, Bs) pump, [Fg ^ b b 

of spiu, - FBlBs\l{‘^f at t/ = 0, (c) + Fi^lBsViW U = L, aud (d) at 

U = oo. The values of the parameters used for these plots are rL=rfl = r, r^ = F^ = 0.1, /3 = O.S/F, uJo ~ fJ- — SF, and 
Br = 0, and all g factors {gl, g^, gg) are —0.44 (bulk GaAs). The hatched areas of (b),(d) denote negative value. 


where iot) = I „. It leads to the BSN curva- 

c \ J d(ixb^) lx=o 


tures 




c=t,4,2 


dll^{a) dpc(a) 
da^ da^ 


— (m -n- n). (74) 


We confirmed the consistency between our results and 
Ref. [s^, which studied the similar system for 0 < U < 
oo using the wideband limit. As we explained at llIIBl 
ipli'Ha) of Eq. dH]) corresponds to — ((/Q(a)|, namely 
In the condition of the wideband limit [i.e., 
Eq. (|6t|l with r(, = 0], we calculated {a) (c =tjii2) 
for 0 < U < oo and confirmed numerically the correspon¬ 
dence between (c =t,4., 2) and — ± ^c'*'(q;)] 

for the charge and spin pump. 

Particularly, in the limit U oo, p 2 vanishes and 
F^^(a) reduces to 


= E 

S=t,i 




da^‘ 


9a" 


where pi°°\a) and are the limits 17 —>■ oo of 

Ps{c() and Is'^(a), respectively. Erom Eq. dZQl) we obtain 










4*.; 




and 


;Moo). ^ ^ r,.,,(l - Mu;,)) - 
" $7 


(76) 

(77) 

(78) 


with LOs = ujQ -\- sgsBs- In the following, we fix a" to 
zero bias (/3h = (3, ph = p) and suppose Eq. (l64ll . Then, 
ls'^(c”)(a) equals v^f{a') given by Eq. (|6T]l and pl°°^(a) 
are given by 


p{sBs) = 




1 + e-7(<^4.-M) _|_ g-P(ojt-F ' 

We emphasize that Fmn°^^ (a) can be obtained by just 


( 79 ) 


(m O njj75) 
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Charge( C/=0) 



Spm(U=0) 


-10 -5 0 5 

/'B-Ds/r 

Charge( U=cc) 



-10 -5 0 5 

(ibBs/T 

Spin(?7=oo) 



-5 0 5 

I‘bBs/^ 


-10 


-5 0 5 


FIG. 2: (Color online)(a) BSN curvature of charge of {Al,Bs) pump, F s^]/^ at U — 0, (b) the BSN curvature 

of spin. [F^l g^ - F^lg^]/^ at [/ = 0. (c) atU^oo, and (d) [F^^g^ - at U = oo. 

The values of the parameters used for these plots are = Bu = F, 7 ^ = F^ = 0.1, and Bl ~ 0 and other conditions are the 
same as Fig. [T] The hatched areas of (a),(c) denote negative value. 



(b) [n'(Bs)-n(-i?s)]-r (c) [n'{Bs) + n'i-Bs)]-r 


0.10 

0.05 

.^0 -20 -10 

/ ■° * 

/ mBs \-0-2 

-30 -20 -10 j 

10 20 30 I' 

\ -qL 



-0.4 

\_^0.10 

-0.5 


f‘uBs 


10 20 30, 


FIG. 3: (Color online) (a) n{Bs) = f{ijJo -b gsBs) (solid line) or p{Bs) (dashed line), (b) n'{Bs) — n'{—Bs), and (c) 


n'{Bs) + n'{—Bs), where n'(±Bs) 


_ J_ dn(B) I 

gs BB 


5 . In all plots, P = 0.5/r, cjq = M and gs = —0.44 x /iB/2. 


replacement. 


Es{a) = f{uJs) -f p{sBs), 


(80) 


in Eq. (|63jl . The charge and spin BSN curvatures of 
{Bl, Bs)t (Al, Bs) pump are given by a replacement 
f'{uJo ± gsBs) —f p'{±Bs) in Eqs. (1551) and (I55)) . where 


p'{Bs) ^'dBs'’ • S™ilar to f/ = 0, the charge and 

spin BSN curvatures of {Bl, Bn), {Al, Ar) pump are 
zero. 

In Figs. [Ila)-I(d), we plot the BSN curvatures of 
{Bl,Bs) pump normalized by (/xb/E)^, where F = E^ = 
Efl and pB = 57.88 /xeV/T is the Bohr magneton. For 
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[/ = 0, the charge and spin BSN curvatures are shown 
in Fig. [Tfa) and Fig. [Hb), and for C/ = oo these 
are shown in Figs. [TJc) and 1(d). The vertical and 
horizontal axes of these plots are the strength of mag¬ 
netic fields Bs, Bl normalized by r//rB- The values of 
the parameters used for these plots are Fl = F/j = F, 
F'^ = F'^ = 0.1, P = 0.5/r, Wo = /i - 3F, Br = 0, and 
g*R = g*R = g*S = —0.44 (bulk GaAs). The BSN curva¬ 
tures of {Ar,Bs) pump normalized by /tb/F are shown 
similarly in Figs. [^a)-2(d). In all plots, = F/j = F, 
7 )^ = = 0 . 1 , Br = 0 , and other conditions are the 

same as in Fig. [TJ In Figs. [T]and[21 the maximum values 
of |r(,(gsi ?5 — gi,Bi,)\/Tb are 0.44 and 0.22 (<1), respec¬ 
tively. The pumped charges and spins are given by Eq. 

(EZl). 

Figure [3Ka) shows the instantaneous average numbers 
of the up spin electron of the QD, n{Bs) = f{oJo + gsBs) 
(for U = 0, solid line) or p[Bs) (for U = 00 , dashed 
line) for /3 = O.S/F, ujq = p — SF, and gs = —0.44 x 
^b/ 2. Because two electrons cannot occupy a QD at 
U = 00 , the magnetic held dependence of p{Bs) is more 
sensitive than f{ujQ+gsBs). Figures (SKb) and 3(c) show 
n'{Bs) ^n'{—Bs) normalized by l/F, where n'{±Bs) = 


In Figs. H a) and 2(c), the charge BSN curvatures of 
{Al,Bs) pump vanish at Bs = 0. This is because the 
hrst term of Eq. (EH) vanishes since n'{Bs) — n'{—Bs) = 
0 for Bs = 0 and the second term vanishes since 
gsBs — gtBb = 0 for Bs = 0 = Br. Similarly, in 
Figs. [IKb) and (d), the spin BSN curvatures of {Br, Bs) 
pump vanish at Br = 0 = Br. The zero lines in these 
plots relate to the cancellation between the hrst and sec¬ 
ond terms of Eq. (I55|) - Figures. [ira),l(c) and Figs. 
[T](b),l(d) are respectively symmetric and antisymmet¬ 
ric under the transformation {Bs,Br) —^ {—Bs,—Br). 
Similarly, Figs. [2Kb),2(d) and Figs. [2{a),2(c) are respec¬ 
tively symmetric and antisymmetric under the transfor¬ 
mation Bs —)> —Bs- We emphasize that pure charge and 
pure spin pumps are respectively realized for {Br,Bs) 
pump and (Aj;,, Br) pump such that the areas S'" in Eq. 
(1671) are symmetric under the above transformations. An 
instance of symmetric area of {Br, Br) pump is a disk of 
which the center is Br = 0 = Br. 

In Wo > /i region, the larger wq — p, the less difference 
between U = 0 and U = 00 becomes. The Coulomb in¬ 
teraction prevents two electrons from occupying the QD. 
This effect is conspicuous in the ujq < p region, although 
it is not important in the ujq > p region. 

As shown in Eigs. [IKa),l(c) and Eigs. [2Ib),2(d), the 
Br dependence of the charge BSN curvature of {Br,Bs) 
pump and the spin BSN curvature of {Ar,Bs) pump at 
U = 0 are more gentle than those at C7 = 00 . It results 
from the behavior of n'{Bs) -|- n'{—Bs) as shown in Eig. 

He). 

As shown in Eigs. [IIb),l(d) and in Figs. [2 K^)j2(c), the 
Br dependence of the spin BSN curvature of {Br,Bs) 
pump and the charge BSN curvature of {Ar,Bs) pump 
are opposite. This is because the leading term (in 


weak magnetic field region) of these are proportional to 
n'{Bs) — n'{—Bs) and its Br dependence is opposite in 
B = 0 and C/ = 00 for wq — /r < 0 as indicated in Fig. 
Hb). This inversion is realized for only wq- /r < 0 region. 
At Wo = p, f'{uJo -f gsBs) - /'(wq - gsBs) vanish. In 
uJo> p region, the signs of f'{uJo+gsBs)-/'{uo-gsBs) 
and p'{Bs) — p'{—Bs) are the same. 

In Figs. [T] and H absolute values of the normalized 
BSN curvatures are smaller than unity. However, we can 
improve this problem by tuning g factors. The first 
and second terms of the right side of Eq. dSH are the 
second and third order in the g factors, and the first and 
second terms of the right side of Eq. (IMll are the first 
and second order in the g factors. If all g factors change 
to —20 (for example for the materials like InAs, InSb), 
the hrst, second, and third order terms become about 
45, 2 000, and 90 000 times. In fact, for these values 
of g factors, the assumption Eq. (IMl) is not appropriate 
for magnetic helds that are not small; we need concrete 
energy dependence of linewidth functions. 


VI. CONCLUSION 

In this paper, we investigated quantum adiabatic 
pump of charge and spin using the FCS-QME (full 
counting statistics with quantum master equation) ap¬ 
proach proposed by Yuge et alM. We studied the 
nonadiabatic effect and showed the correspondence be¬ 
tween our approach and the real-time diagrammatic 
approach^f ^ill B|) . and conhrmed the consistency be¬ 
tween the two methods in the concrete model, the one 
level interacting quantum dot (QD) ( giV Cl and ilVCl) . 
Moreover, in ijll B[ we showed that the Berry-Sinitsyn- 
Nemenman (BSN) phase derived under the “adiabatic” 
condition (which makes the Berry phase like treatment 
appropriate) has the nonadiabatic information. The 
FCS-QME picks out one higher order nonadiabatic piece 
of information from the solution of the QME. Partic¬ 
ularly, the instantaneous steady state (the zeroth or¬ 
der of the pumping frequency) gives first order response 
(pumped current). This fact may be related to Ref. [1^ 
which is connected the BSN vector and the liner response 
theory of the QME. 

We generalized the ECS-QME approach to the multi¬ 
counting field f ijll AD and studied the QDs system weakly 
coupled to leads {L and R) in ^IVI and (|Vl We showed 
that the pumped charge and spin coming from the instan¬ 
taneous steady current are not negligible when the ther¬ 
modynamic parameters (the chemical potentials and the 
temperatures of leads) are not fixed to zero bias in nonin¬ 
teracting QDs fi lIV Bp and an interacting QD f ilV BD . To 
observe the spin effects, we consider collinear magnetic 
fields, which relate to spins through the Zeeman effect, 
with different amplitudes applying to the QDs {Br) and 
the leads {Br and Br). We focused on the dynamic pa¬ 
rameters {Br, Br/r and the coupling strength between 
QDs and leads, Ar/r) as control parameters. 
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In one level QD with the Coulomb interaction U, we 
analytically calculated the BSN curvatures of spin and 
charge of {Bl,Bs) pump and {Al,Bs) pump for the 
noninteracting limit (U = 0) and the strong interaction 
limit (U = oo) using the rotating wave approximation 
(RWA) defined as the long coarse-graining time limit of 
the coarse-graining approximation (CGA). The difference 
between U = 0 and U = oo appeared through the in¬ 
stantaneous averages of the numbers of the electron with 
up and down spin in QD. For {Bi^,Bs) pump, the en¬ 
ergy dependences of linewidth functions, which are usu¬ 
ally neglected, are essential. Additionally, the adiabatic 
modulations of (Bl^Bh) or (A^, A/j) can pump neither 
charge nor spin. 

In this paper, only U = 0 and U = oo limits are stud¬ 
ied. In fact, we can analyze finite U based on Eq. dZH). 
This is our future work. Recently, Yoshii and Hayakawai^ 
studied adiabatic pump of charge by only the thermody¬ 
namic parameters using the same approach in a similar 
system without magnetic fields (for finite U). The work 
(for the thermodynamic parameters) involving the above 
problem is contrastive to our work (for the dynamic pa¬ 
rameters). 
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a 
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a 
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Hnm,kl — {Hs)nkSl m ^nk (All) 

Bnm,kl ~ '^^{^a)nk{Ba)lm- (^ 12 ) 
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Appendix A: Liouville space 


By following correspondence, an arbitrary linear op¬ 
erator (which operates to the Hilbert space) • = 

• \'m)\n)(m\ is mapped to a vector of the Li¬ 


ouville space^i^, !•)) = 

\n){m\ 

Tr{\m)(n\n'){m'\) 

Tr(Al'H) 

Tr(.) 


,min\ • \m)\nm)): 


i —|nm)). 

(Al) 

^—> ({nm\n'm')), 

(A2) 

^ ((A|H)), 

(A3) 

^ ((1|.)). 

(A4) 


Here, {|n.)} is an arbitrarily complete orthonormal ba¬ 
sis. The inner product of the Liouville space is defined 
by the Hilbert-Schmidt product [Eq. (IA3I) ]. The Her- 
mitian conjugate of !•)) is defined as ((•! = (l*)))^ = 
Snm(^l • |w)*((nm|. An arbitrary linear superopera¬ 
tor J which operates to any operator (•) is mapped to a 
corresponding operator of the Liouville space (J) as 

I J-)) = J|-)). (A5) 

The matrix representation of J (or J) is defined by 

Jnm,kl = ((nml J|fc^)). (A6) 


Einally, we consider the current operators defined by 
Eq. kf^ = is given by 

= (A13) 

a 

Hence the current operators defined by Eq. (HD) are given 

by 

(A14) 

a 


Appendix B: Details of FCS-QME 

is governed by the modified von Neu¬ 
mann equation^: 

= (Bi) 

where Htot{t) is the total Hamiltonian and [A,= 
A^B - BA_^ with A^ = 

is given by = Hsit) + Y,b{^b{t) + Hsb{t)], 

where Hs is the system (denoted by S) Hamiltonian, 
Hb is the Hamiltonian of bath b, and Hsb is the tun¬ 
neling interaction Hamiltonian between S and bath b. 
In the following, we suppose Eqs. (Pl) . dMl) and dUl) 
for an arbitrary number of leads (b = 1,2, ••• ,M). 
The initial condition of is given by^ Ptot(O) = 
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-P{o^}Aot(0)P{o„}- Here, {o,^} denotes eigenvalues 
of {O^} and P{o^} is a projection operator defined by 
P{o,^}P{oi} = ^{ou}Y[i_i^oij,,o'^^ and 
suppose ptot(O) = p(0) 0 psiao) 

where 

PbK) = <S) (B2) 

b 

with Sf,(at) = Trf,[e“^'’^*dBi,(t)-M6(b^i.]] a^d Trj, de¬ 
notes the trace of lead b. Then, Ptot(O) = P(0) ® 
PMPB{ao)P{o,} obeys. If all are given by 

T.b,k,crObka4kaCbka with real numbers P^o^y com¬ 

mutes to p_b( 0) and Ptot(d) = p{0)^Pb{o:o) obeys because 
S{o^} = ^■ 


Now we move to the interaction picture. An oper¬ 
ator in the interaction picture corresponding to A(t) 
is defined by (t) = t/J(<)A(t)f/o(i) with = 

-i[Hs{t) + Y.bHb{t)]Uo{t) and t/o(0) = 1. The sys¬ 
tem reduced density operator in the interaction picture 
is given by p^’^{t) = Trieads[Pto?(0] where = 

Uo{t)pl, {t)UQ{t) and Trieads denotes the trace of lead 
space, ploti^) is governed by 

^PlM = .,[HUt).pioUt)k, (B3) 

with Hint = ^b^Sb- Up to the second order perturba¬ 
tion in Hint, we obtain 


/ t+rcG 

du I ds Trieads [^4 («): WPS («*)]%] X } 


= -h TcgU^cg 


using the large-reservoir approximation pi^{t) « 
p^’^{t)® pB{at) and Trieads(a*)] = 0- The ar¬ 
bitrary parameter tcg (> 0) is called the coarse-graining 
time. The coarse-graining approximation^Si^ (CGA) is 
defined by 

lpP^t)=L^^^{t)p^’Ht). (B5) 

If the cycle time of the modulation of control parameters, 
T, is much longer than the coarse-graining time tcg, the 
superoperator L^cg i® described as a function of the set 
of control parameters at time t. In this paper, we suppose 


(B4) 


T ^ TCG- Moreover, tcg should be much shorter than 
the relaxation time of the system, ts ~ i. On the other 
hand, Tg <C t should hold for the adiabatic pump. Hence 
Tcg "C ^ <C t should hold. 

In the Schrodinger picture, Eq. (IB5I1 is described as 

= -^[Hs{t),p^t)]+H>^^^{at)p>^{t). (B6) 

Here, the superoperator operates to an operator 

• as 


n?cGH* = m] II 


(X, tcg, W, (w') • [a„s(w)]t - (tcG, W, w ') • (w)] (w') 


oj.oj' n,m,s,s' 


-\‘^b,ns,ms' (tcg , W, w') [ans (w)]^ams' (w') • (X, TCG, W, Uj') [Ums' (w')]^ * ^ns (w) 

-^^ftns.ms'(^CG,w,a;') •ans(w)[a„s'(w')]^ - i$+„^^^^,(TCG,W,w')ans(w)[ams'(w')]^ • 

H [-\'^b,n,,ma'i^CG,UJ,Uj')[ans{uj)]'<ams'{uj') + ^ (tcG , W, ^Oons (w) [Oms' (w')] ^'IBT) 


uj.uj n,m,s,s 


where 


A^(x,TcG,cd,cd') = / ^»-^^(^^x)TCGsinc( ^^^^^ ‘^^ )sinc( "'^°^^^ (B8) 


and A=‘=(tcg,w, w') = A=‘=(x = 0, tcg, w, w'). Here, sinc(a;) = sinx/a; and A=‘=(H,x) denotes one of (H, x). 













15 


'^b,ns,ms'i^^x), where 


k,a 

^b)ns,ms'(^’= 27r V5,j,^ „gVbfeo-,ms'/b(ebfccr)e e S^Ebka — XI), 

k,a 

'^b,ns,mA^^X) = 2Y,Vbka,nsVbA^^A^ ” /b (^bfc. )]e*^-P - 

k,(T 

Kns,ms'iAx) = 2^P,L,„,Pbfc<.,ms'/b(ebfc<.)e-*>^-e-*^-^-P- 
k,(T 

I- 


1 

^bk<7 XI 
1 

^bka XI 


(B9) 

(BIO) 

(Bll) 

(B12) 


with Vbka,ns = A^Vbk(T,ns- Xba and \ba denote the 
counting fields for Nbn and ^k^bkc\f.^Cbka respectively. 
/b(e) = [exp(/3{,(/r{, — e)) + 1]“^ is the Fermi distribution 
function and P denotes the Cauchy principal value. The 
eigenvectors UnsA) are given by 

ansA) = ^^u:fic.,i^AA{Ea\ans\E0){EA (^13) 

Oi,j3 

with ui/sa = Ep — Ea and E[s\Ea) = Ea\Ea)- bo is one of 
the elements of {u}pa\ {Ea\ans\Ep) ^ 0 ^(n, s)}. 

We explain the rotating wave approximation (RWA). 
If Hs is time-independent, factors )‘ appear in 

Eq. (IB16F The usual RWA^^ is approximating 
as in Eq. (IB16F However, if Els is time de¬ 

pendent, the generalization of the RWA is unclear. In 
this paper, the RWA is defined as the limit tqg oo 
of the CGA. In this limit [tcg ' inintj^^/(|w — w'l) ^ 
1], X^{x,tcg,ui,uj') Ki XAb^x)A,u,' holds because of 
the fact that limT-^^-^oo TcGsinc ^'^°^^~'^^ sinc ^'^'^^^~‘^ ^ = 
2'K5^^i^i8{Xl — w)- If Hs is time independent, this RWA is 
equivalent to usual RWA. 

In the Born-Markov approximation (without RWA), 
Eq. (IB16I) sometimes violates the non-negativity of 
p[t){= p>b=°{t))^. The QME of the RWA or the CGA 
is the Lindblad type (see Eq. (IB7I1 at x = 0) which 
guarantees the non-negativity^. In QME of the RWA, 
the diagonal and off-diagonal elements of p{t) are de¬ 
coupled. Sometimes, the decoupling is justified by the 
superselection rules^^— : if two states in the system 
I to) and |n) differ in an observable which is conserved 
in the total system, unavoidable interactions lead to 
rapid decay of (TO|p(t)|n). The standard example is 
the electron number which is conserved for QDs cou¬ 
pled to non-superconducting leads. The total spin pro¬ 
jection is also conserved for unpolarized or collinearly 
polarized leads when spins do not flip by the interac¬ 
tion Hsb- Hence, in our setting, Eqs. (PI). (P|) with 
Vbka,ns oc Sa,s, and Ep. lHBl) . which was used after (II V Cl 
the decoupling between the diagonal and off-diagonal 
part of (TO.|'TO 4 ,|p(t)|n|n 4 ,) is justified. However, in the 
general system (discussed until (IIVBI) . the decoupling 
can not be justified by the superselection rules. In fact, 


as explained the following, the CGA can break the de¬ 
coupling. However, the QME without decoupling is too 
hard to analyze analytically. 

The FCS-QME of the CGA (for finite tcg) is more 
difficult to analyze than that of the RWA. For the one 
level system Eq. (H51) . the matrix representation of the 
Liouvillian of the RWA by the number states {|n.|-n^)} 
is block diagonalized to the diagonal part (spanned 
by {|ntn 4 ,)(n|n 4 ,|}„.f,„^=o,i) and the off-diagonal part 
(spanned by {\n-fni){m^mi\}(n^n^{rnt,rnA- The di¬ 
agonal block is given by Eq. (|65)) and the off-diagonal 
block is a (12 x I2)-diagonal matrix. However, one of the 
CGA has off-diagonal components which permit transi¬ 
tions between diagonal and off-diagonal if spins can flip 
by tunneling Hint- This is consistent with the above dis¬ 
cussion about the superselection rules. Even if spins can¬ 
not flip, off-diagonal block is not diagonal. Particularly 
for U = 0, the Liouvillian dose not reduce to a summa¬ 
tion of one particle Liouvillian [Rr^(a)0i4,-|-l'|'0lF^(a)]. 
The study of differences between the RWA and the CGA 
is a future work. 

Finally, we recognize the Born-Markov approximation 
(without RWA). From Eq. (IB3I) . we obtain pl'AA = 

Pto?(0) -Ao^u [H-AA), pAA)]x- Substituting this to 
Eq. (IB3p . we obtain 


^dM = .i[HAit),pA{o)]x 

-WLit), f du [HA{u),pIAA)]x]x- (B14) 

Moreover, substituting Ptoti'^) = Ptot(0 “ 
i^ds [HAis),plAis)]x to Eq. (|BI4| and using 

pIAA ~ P^’^ApsAt) and Trieads[i?4(OPB(ao)] = 0, 
we obtain 


‘''■'■’‘W = - /b. Tr„. 


dt 


{[HLA^ [hLA), P^’^ApbAAxIx}^ (B15) 


up to the second order perturbation in H-mt ■ Because the 
integrand decay as .^^here tb is the relaxation 
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time of the baths, we can replace du to J^^du, and 
we obtain 

dp^’^{t) 


dt 


ds Trieads 


{[HL(t), [HLt(t - s),/’^(i)ps(at)]x]x}- (B 16 ) 


This is called the Born-Markov approximation (without 
RWA). Equations (IB15I) or (IB16I) are also derived from 
Eq. (IB4I) . In the second term of the right side in the 
first line of Eq. (lB4l) . we can replace ds to ds or 
f^^ds. After this replacement, by taking arcc Ir ,-o’ 
we obtain Eqs. (|B15p or (jB16l) respectively. In the Eqs. 
(IB15P or (jB16() (after moving to the Schrodinger picture), 
we can replace {au}o<u<t or {at_s}o<s<oo to at because 
of tb t. 


Appendix C: The time evolutions of c^{t) 

In this section, we derive the time evolution equations 
of c^(t) of Eq. (0. The left hand side of the FCS-QME, 
ilP^(t))) = K^{at)\p^{t))), is 


I dt 

n 

+c)((t)e^"(‘)A„(aO|p(((at))) 

+c)((t)e^"W||p^(a,)))}- (Cl) 


And the right hand side of the FCS-QME is 

k^atWm =Y.cme^^^^^'>k^at)\pl{at))) 

n 

= ^c^(t)e^^(‘)A„(aQ|p)((aQ)).(C2) 

n 

Hence we obtain 


=0. (C3) 

, and using 


Applying {{l^[at)\ to Eq. 
((^n(«)|Pm(«))) = ^6 obtain 




dt 


By the way, the time derivative of Eq. ©, 
k^{at)\p^{at))) = AX(aQIpX(at))), is 

dk^{at)^ d\pl{at))) 

+ K^a,) ^ 

= + (C5) 


Applying {{kn{^t)\ to this equation, we obtain 

nix MI^x_L \x \inx \i (o-^))) 
\\^m(o^i)l \Pn\^t)l} A Ay,.j (at)((/^(at)I 

_ I \x (^ M^l/^n(^i))) f 

Omn + •^n('^t)((^m('^t)| j (C6) 


dt 

and it leads to 




d\pUat))) {mat)\^^^^\pUat))) 


dt 


Am (at) - Xn{at) 


{,C7) 


for AX (at) ^ X^{at). Substituting this to Eq. (IC4E we 
obtain 

^^=-((lx(at)|||pX(at)))cX(t) 


+ E 

n(^m) 


AS(t)-A^(t)((^^(«*)l^^|P^(«t)))t 


Am(at) - Xn{at) 
The above equation can also be written as 


(C8) 


dt 


= E 




n{^m) 


(C9) 


{mat)\^i^^\pUat))) 

"" XUat) - AO(at) ’ 

where c^(t) = with 

Vmit) = {{1■^{as)\■^J\Pl^{as))) 

= (CIO) 

Here, C is the trajectory from ap to at, a^ are the fcth 
component of the control parameters, and p^{t) = 0{1) 
since ((^x (at)|^|pX (at))) = 0(uj) with w = 27r/r. In 
the right hand side of Eq. (ICQl) . the dominant term is 
n = 0 if m 0 because ReAg(a) > ReAx(a). Using 

= o(ra;), Ax(at) = 0(r), = 0(1} and 

Cg(t)e^o = 0(1), we obtain 




dt 




(Cll) 


and 


c^{t)e^^d) = o(^ r ds 
Jo 

For X = 0, Eq. (IC12I) is also derived from 


(C12) 


/o“(0 = PiO - Po(at) = E Cm(Qe^'"(*Vm(at):(C13) 

771^0 

and Eqs. (USD and (EZD. 
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Appendix D: The validity of the adiabatic expansion and rimax = [^max] with 


In the derivation of the QME with CGA, when going 
from Eq. (IB5I) to Eq. (IB 6 I) . we used the following type 
of approximation: 



ds s, u; t) 



/ U 

ds G'([at]; s, u; t). 


(Dl) 


Here, G([a]“, s, it, t) ^ [q,jm _ ( 04 /)“,^^ 

is the control parameters trajectory and [a*] is the tra¬ 
jectory which at' = at {s < t' < u). Similarly, in the 
Bron-Markov approximation (BM), when going from Eq. 
(IB15I) to Eq. dll), we used 


du G{[a\l;u,t;t) 


du G([at];It, t; t). (D2) 


Considering the corrections of the above approximations, 
the QME are given by 

^^=/C(t)|p(t))), (D3) 

/C(t) = k{at) + k\t), k'{t) = 0{VuJTx), (D4) 


with Lo = 27r/T and tx = tcg for CGA; tx = tb for 
BM. The corrections are also discussed in Ref. [3l|. The 
discussions between Eq. (l24t and Eq. 1221 ) are cor¬ 
rect after replacing K{at) —>■ /C(t), TZ{at) —>■ 'R-{t) and 
Po{oit) —t Po{t)- Here, po{t) and 'R-[t) are defined by 
K.{t)\pt){t))) = 0 and 'R-{t)lC{t) = 1 - |/5o(t)))((l|, respec¬ 
tively. Equation (1251) is corrected to 




dt 


dt 


00 , 00 

E [^W|]"|poW)) ^ E IP“^”Hi)))(D5) 


with p“(t) p{t) — po{t)- The corrections are given by 
po = po[l + 0{^Tx)l n = n[\ + 0{^Tx)l (D6) 

and 

/5“(")W -P“(")W = 0((^Ywtx). (D7) 

Next, we consider the reasonable range of n of 
Because = G(^)" and poit) — po{at) = 0 {ujtx), 

the reasonable range is n < where n^ax is deter¬ 

mined by 

(p) < wry < (-) . (D8) 

Let us consider that reasonable concrete values of the 
parameters in our model f llHIl) : oj = 10^ MHz, T = 
10 peV=0.116 K, 1/r = 65.8 ps, tcg = 1 ps, and 
Tb =0.1 ps. These values lead to 

curcG = 10-®+^ cutb = 10-"+P, ^ = 10-418+p,(D9) 


«=GA). (BM). (DIO) 

Here, [n] means the biggest integer below n. At p = 0, 
fimax = 1-44 (CGA), 1.67 (BM) and at p = 3, nmax = 
2.54 (CGA), 3.39 (BM). The larger the nonadiabaticity 
(p), the larger Umax becomes. 


Appendix E: Proof of Eq. (13111 

First, using Eqs. (l24l) and (l23l) . we obtain 
{{l\W^{a)TZ{a)k{a) = ((l|IE^(a) - A^(a)((l|. Next, 
((1 o(q()| = ((1|, Ao(a) = 0, and Eqs. Q and (l 22 l) 
lead to ((/(((a) I AT (a) = Ag(a)((l| — ((l|H2(a). Hence 
[((l|lT^(Q;)77.(a) -I- ((1q (q;)|] Ar(Q;) = 0 and it leads to Eq. 

To prove Eq. (1^ only Eq. (IMl) is required and 
k{a)TZ{a) = 1 — |po(Q()))(( 1| is not necessary. Addition¬ 
ally, the pseudoinverse of the real-time diagrammatic ap¬ 
proach Eq. (15^ satisfies 

E - P^ ^ E (El) 

K K 

which corresponds to our 

n{a)k{a) = 1 - |po(a)))((l| 7 ^ k{a)n{a). (E2) 

Equation (IMI) is shown also as follows. Eq. (I24|) and 
((l|A'(a) = 0 lead to k{a)TZ{a)k{a) = k{a), which 
implies 

A(a)7^(a) = l-|a(a)))((l|, ((l|a(a))) = 1. (E3) 

Applying ((1| to Eq. (IMl) . we obtain ((l|77.(Q;)Ar(a) = 0, 
which is equivalent to 

((l|7^(a)=C(a)((l|. (E4) 

By the way, differentiating Eq. o for n = 0 by iXfi, we 
obtain 


{{li;{a)\k{a) + {{l\k>^{a) = ((l|A(;(a). (E5) 

Applying TZ{a) to this equation and using Eqs. (IE3I) and 
(IE4I) . we obtain^i^ 

{{li;{a)\ = -{{l\k^^{a)n{a) + c^(a)((l|, (E 6 ) 

c,{a) = C(a)A(;(a) + C(a)|a(a))). (E7) 

Eq. (IE 6 I) becomes Eq. (I^Tl) because of Eq. (1^ . Partic¬ 
ularly, Yuge^ used 

7^(a) = -lim rdte^fo)*(l-|po(a)))((l|), (E 8 ) 

which satisfies Eqs. (IE3p and (|E4I) with cr{a) = po(a), 
C(a) = 0 and Eq. (IMl) (in Ref. [Sg, C(a) was incorrectly 
set to — 1 ). 
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